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INTRODUCTION 


This memo documents the linear aerodynamic and thrust perturbation model of the 
X-29 simulation model resident in the Control Integration and Assessment Branch of 
Wright Laboratory (WL/FIGD) at Wright-Patterson Air Force Base, Ohio. This model 
was used to perform flying qualities evaluations of the X-29 aircraft at Wright 
Laboratory during the airplane’s development in the early to mid-1980s. 


This work was conducted in order to document the X-29 longitudinal and lateral- 
directional linear aero model at numerous trim flight conditions for use in future 
research work in flying qualities. The X-29 is representative of typical fighter aircraft, 
even with its forward-swept wing, although the airframe is very unstable longitudinally. 


This memo is divided into four sections: Introduction, Definition of Terms, Explanation 
of Computer Printouts, and Computer Printouts sections. Definitions of all terms 
associated with the linear model are given in tie Definition of Terms section. 
Explanation of all of the data in the computer printouts is given in the Explanation of 
Computer Printouts section. The Computer Printouts section contains the actual 
computer printouts for each trim flight condition. 


Mr Dan Goddard of WL/FIGD wrote much of the original computer software that 
generated the linear model given in this memo. The software was modified where 
necessary to make the computer printout easier to understand and to complete the 
lateral-directional section of the software. Mr Brian Stadler and Mr Joe Nalepka of 
WL/FIGD resurrected and documented the X-29 aero model for this project and 
helped immensely in explaining the computer system and software at WU/FIGD. 


DEFINITION OF TERMS 


The aerodynamic forces and moments for the model are defined in terms of 
nondimensional aero coefficients, C,, C,, C,, and C,, C,,, C,, respectively. The x axis 
is positive out the front of the airplane, the y axis is positive out the right wing, and the 
z axis is positive out the bottom of the airplane. The aerodynamic forces and 
moments are functions of twelve variables, seven in the longitudinal axis and five in 
the lateral-directional axis. These twelve variables are: angle of attack, pitch rate, 
Mach number, angle of attack rate, strake deflection, canard deflection, flaperon 
deflection, angle of sideslip, roll rate, yaw rate, aileron deflection, and rudder 
deflection. Thrust variation is also included in the linear model and is a function of two 
variables, Mach number and altitude. Thrust is aligned along the aircraft centerline 
parallel to the x axis (no inclination angle) above the center of gravity (negative 
pitching moment). 


The linear model is produced by perturbing each of the twelve variables, one at a 
time, in a positive direction and then in an equal negative direction from the trim value. 
The partial derivatives required for the nondimensional and dimensional derivatives are 
calculated by dividing the total variation in each nondimensional aero coefficient or 
thrust by the total variation in each of the twelve independent variables. Altitude is not 
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varied in these linear models. All angular and angular rate independent variables are 
varied plus and minus one degree or one degree/second, respectively. Mach number 
is varied plus and minus 0.05. 


All aero coefficients and nondimensional and dimensional derivatives follow the 
definitions given in Appendix A of NASA CR-2144 (Reference 1) except for two 
variations. It was decided to make all sign conventions agree for forces and 
moments. This is so all nondimensional coefficients and nondimensional and 
dimensional derivatives are positive in the positive directions of the axis system. 
Therefore, C, is the opposite sign of the C, in NASA CR-2144 and C, is the negative 
of the C,, in NASA CR-2144. The lateral-directional L and N dimensional derivatives 
aré converted to primed derivatives as defined in NASA CR-2144 for print out. Please 
note that all coefficients and derivatives are in a body axis system, not a stability axis 
system. 


The dimensional derivatives are used to generate a longitudinal and lateral-directional 


state space model according to the following equations. The general form of the state 
space model is: 


X=Ax+Bu 


where A is a 4x4 matrix and B is a 4x3 matrix for the longitudinal model or a 4x2 
matrix for the lateral-directional model. 


Longitudinal State Space Model 


The state vector x contains forward speed, angle of attack, pitch rate and pitch angle. 
The input vector u contains strake, flaperon and canard deflections. The elements of 
the A matrix are: 


MZ. x,Z, 
Ss +. anu A 1,2 : i 
A(1,3)=Xq- We ere A(1,4)=-9 cos tyr iF > 
0° a 0-s 
2, 2, 
= A(2,2)= 
A(2,1) U,-Z, (2,2) Wz, 


q 


02a 


A(2,3)= me 


M,Z; 
Uy-Z, 


A(3,1)=M* + 


A(3,3)=Mg+ 


A(4,1)=0 


A(4,3)=1 


The elements of the B matrix are: 


x,Z 
B(1,1)=X_ + 
- U,-2,, 


M;Z5, 


B(3,1 )=My + U,-Z, 


X25, 


Uy -Ly 


B(1,2)=Xq+ 


M.Z, 
B(3,2)=M, + aA, 
= 


UZ, 


M,( Up +Z)) 
Us-Z, 


-g siN®, 


A(2,4)=—-_-- 
(2.4) 


M,Z, 
Uy-2, 


A(3,2)=M,+ 


A(4,2)=0 


A(4,4)=0 


Zs, 


B2,1)=— ¥ 
poe, 


B(4,1)=0 


2, 


Ty-Z 


B(2,2)= 


B(4,2)=0 


X, 
B(1,3)=X,,+ — 


B(2,3)= 


- 


B(3,3)=M,, + 


M.Z, 
= B(4,3)=0 
= 


U,-Z, 


a 


Terms with subscripts of 5,, 5, or 5, are associated with the strake, flaperon, or 
canard control surfaces, respectively. 


Lateral-Directional State Space Model 


The state vector x contains angle of sideslip, roll rate, yaw rate, and bank angle. The 


input vector u contains aileron and rudder deflections. The elements of the A matrix 
are: 


A(i tege A(1,2)= VS 
M19 A(t 4a 
A(2,1)=Ly A(2,2)=L) 
A(2,3)=L) A(2,4)=0 
A(3,1)=N5 A(3,2)=Ny 


A(3,3)=N; A(3,4)=0 


A(4,1)=0 A(4,2)=1 
siné 

A(4,3)= 2 = 

(4,3) <56, A(4,4)=0 


The elements of the B matrix are: 


B(1,1)=-% BU1,2)=-% 
Vz. V7, 
B(2,1)=Ls' B(2,2)=Ly 
B(3,1)=Ng_’ B(3,2)=N,’ 
B(4,1)=0 B(4,2)=0 


Terms with subscripts of 5, or 5, are associated with the aileron or rudder control 
surfaces, respectively. 


EXPLANATION OF COMPUTER PRINTOUTS 


The computer printout is divided into six pages for each trim condition. The first page 
is a printout of the trim condition including the aircraft state, positions of the controls, 
the trim aerodynamics, and aircraft constants. The second page is a printout of all of 
the nondimensional aero derivatives and thrust-related derivatives. The third page is a 
printout of the longitudinal dimensional stability derivatives, the longitudinal state space 
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model, the roots of the longitudinal characteristic dynamic equation, and the period, 
damping ratio, damped and undamped natural frequencies for each complex pair of 
roots present in the characteristic equation. The fourth page is a printout of all 
longitudinal transfer functions for each state output to each control surface input and 
the transfer function for normal acceleration output at the center of gravity to each 
control surface input. The fifth page is a printout of the lateral-directional dimensional 
derivatives, the lateral-directional state space model, the roots of the lateral-directional 
characteristic dynamic equation, and the period, damping ratio, damped and 
undamped natural frequencies for each complex pair of roots present in the 
characteristic equation. The sixth page is a printout of all lateral-directional transfer 
functions for each <tate output to each control surface input and the transfer function 
for side acceleration output at the center of gravity to each control surface input. 


The first page summary of the trim condition is fairly self-explanatory. All trim 
conditions are for straight, wings-level, 1g flight. All angular values are in degrees. 
The variable "VTRFPS” is the aircraft trim airspeed in feet per second. The fuselage 
station and waterline position of the center of gravity are given as the variables "XCG" 
and "ZCG", respectively, and are in inches. The flight control system MODE status 
logicals at the bottom of the AIRCRAFT STATE section are a carryover from the 
original linearization software code and have no effect on the linear model given in the 
printout. The model is strictly for the airframe aerodynamics and thrust variation with 
Mach. Based on the signs of the nondimensional and dimensional derivatives, the 
three longitudinal controls, canards, flaps, and strakes are positive trailing edge down. 
The rudder is positive trailing edge left and the ailerons are positive right wing surface 
trailing edge up. The thrust at trim is shown in pounds. Also, although not printed 
out, the thrust line is above the center of gravity by a distance of 75.204 inches minus 
the value of "ZCG" (negative pitching moment). The variable "% RPM" is the percent 
of full throttle required to set the thrust at the trim value. The nondimensional aero 
coefficients at trim are given in the DYNAMICS section and are positive in the positive 
directions of the axis system as explained earlier. 


The second page of the printout is a listing of all of the nondimensional aero 
derivatives and thrust variation with Mach perturbation. The calculation of the 
nondimensional derivatives and thrust variation is consistent with Appendix A of NASA 
CR-2144. 


The third page of the printout lists the longitudinal dimensional stability derivatives, 
with units shown. The thrust variation with Mach perturbation has been included in 
the calculation of these dimensional derivatives. The longitudinal state space model A 
and B matrices are then printed out. All angular states, "ALPHA", "Q", and "THETA", 
and the control surface inputs, "STRAKE", "FLAPRN", and “CANARD? in this model 
have units of radians or radians/second, as appropriate, and the speed state, "VEL", 
has units of feet/second. The system roots of the A matrix are then printed out in 
realimaginary form. Finally, the period, damping ratio, and damped and undamped 
natural frequencies for any complex pair of system roots is printed out. 


The fourth page of the printout lists the transfer functions in Laplace form for the 
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output of each state to each control surface input. The units on the transfer functions 
are the same as the units on the states and surface deflections in the state space 
model. Also, the transfer function of the acceleration at the center of gravity in the 

z direction to each control surface input is also listed as "AN/STRK", "AN/FLAP", AND 
"AN/CAN". These transfer functions are based on the states ak to the 
following equation: 


AN=U,6-Upq+(gsin®,)@ 


The units on this acceleration are feet/second? and it is positive in the positive z 
direction. 


The fifth and sixth pages repeat the same type of information for the lateral-directional 
motion as the third and fourth pages did for the longitudinal motion. The lateral- 
directional dimensional derivatives have units as shown. The lateral-directional state 
space model states and surface deflections all have units of radians or 
radians/second, as appropriate. The transfer functions include the acceleration at the 
center of gravity in the y direction for aileron and rudder deflections, "AY/AIL” and 
"AY/RUD", respectively. The transfer functions are based on the states according to 
the following equation: 


AY=UpB-(g coS Oo) + Upr- Wop 
The units on this acceleration are feet/second’ and it is positive in the positive 
y direction. 
Computer Printout Variables By Page 
Page 1 


The variables printed out on page 1 are as follows: 


ACGW Aircraft gross weight, pounds 

AREA Reference wing area, feet” 

AILERON Trim aileron deflection, positive right wing surface trailing edge up, 
degrees 

ALPHA Trim angle of attack, positive nose up, degrees 

ALT Altitude above sea level, positive up, feet 


AR MODE 
BETA 
CANARDS 
CBAR 
CLMAERO 


CMMAERO 


CNMAERO 


CXAERO 
CYAERO 


CZAERO 


DR MODE 
DYN PR 
FLAPS 
GAMMA 
GEAR DOWN 


NORMAL MODE 


Flight control system mode, not used 

Trim angle of sideslip, positive nose right, degrees 

Trim canard deflection, positive trailing edge down, degrees 
Reference chord, feet 


Trim nondimensional x-axis aero moment coefficient, positive right 
wing down 


Trim nondimensional y-axis aero moment coefficient, positive nose 
up 


Trim nondimensional z-axis aero moment coefficient, positive nose 
right 


Trim nondimensional x-axis aero force coefficient, positive forward 


Trim nondimensional y-axis aero force coefficient, positive out the 
right wing 


Trim nondimensional z-axis aero force coefficient, positive out the 
bottom of the airplane 


Flight control system mode, not used 

Trim dynamic pressure, pounds/feet” 

Trim flaperon deflection, positive trailing edge down, degrees 
Trim flight path angle, positive up, degrees 

Landing gear up if F(alse), down if T(rue), always F 

Aircraft x-axis moment of inertia, slug-feet® 

Aircraft product of inertia, slug-feet* 

Aircraft y-axis moment of inertia, slug-feet” 

Aircraft z-axis moment of inertia, slug-feet” 

Mach number, dimensionless 


Flight control system mode, not used 


NX 


NY 


RUDDER 
SPAN 
STRAKES 
THETA 
THRUST 
VC KTS 
VTRFPS 
XCG 


ZCG 


Trim forward acceleration in g’s, positive forward, dimensionless, 
always zero 


Trim side acceleration in g's, positive to the right, dimensionless, 
always zero 


Trim normal acceleration in g's, positive out top of airplane, 
dimensionless 


Trim throttle setting, percent of full setting 

Air density, slugs/feet® 

Trim rudder deflection, positive trailing edge left, degrees 
Reference wing span, feet 

Trim strake deflection, positive trailing edge down, degrees 
Trim pitch attitude, positive nose up, degrees 

Trim engine thrust, positive forward, pounds 

Trim calibrated airspeed, knots 

True airspeed, positive forward, feet/second 


Fuselage station position of the center of gravity, positive forward, 
inches 


Waterline position of the center of gravity, positive up, inches 


Page 2 


This page is fairly self-explanatory. All nondimensional stability derivatives are per 
radian for angular perturbations and per Mach number for the Mach derivatives. Thrust 
is separated from the aerodynamic derivatives as shown in the far right column. 
Thrust only varies with Mach number. Thrust variation with all other longitudinal 
perturbations is zero as shown. Thrust variation with Mach number has units of 
pounds per Mach number as shown. 


Page 3 


All dimensional stability derivatives have units as shown at the top of each column and 
at the left of each row. For example, the units for Z, are feet/second” per radian. The 


"feet/second”” 


is shown at the top of the "Z” column just under the "Z” and the "per 


radian” is shown at the left next to the "ALPHA". 
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The state space model variables are as follows: 


ALPHA 
CANARD 
FLAPRN 


VEL 


Angle of attack, positive nose up, radians 

Canard deflection, positive trailing edge down, radians 
Flap deflection, positive trailing edge down, radians 
Pitch rate, positive nose up, radians/second 

Strake deflection, positive trailing edge down, radians 
Pitch attitude, positive nose up, radians 


X-axis velocity, positive forward, feet/second 


The SYSTEM POLES are the roots of the characteristic equation of the A matrix in 
real, imaginary form. 


The DYNAMICS variables for each pair of complex SYSTEM POLES are as follows: 


PERIOD 
WD 

WN 
ZETA 


Period for one cycle, seconds 

Damped natural frequency, radians/second 
Undamped natural frequency, radians/second 
Damping ratio, dimensionless 


Page 4 


This page lists all of the longitudinal transfer functions in polynomial form for each 
state output to each contro! surface deflection input. It also gives the transfer function 
for the z-axis acceleration at the center of gravity to each control surface input. The 
output and input variables are as follows: 


ALP 
AN 


Angle of attack, positive nose up, radians 


Z-axis acceleration at center of gravity, positive out bottom of airplane, 
feet/secon 


Canard deflection, positive trailing edge down, radians 
Flap deflection, positive trailing edge down, radians 


Pitch rate, positive nose up, radians/second 
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STRK Strake deflection, positive trailing edge down, radians 


THT Pitch attitude, positive nose up, radians 
U X-axis velocity, positive forward, feet/second 
Page 5 


All dimensional stability derivatives have units as shown at the top of each column and 
at the left of each row the same as page 3. 


The state space model variables are as follows: 


AlL Aileron deflection, positive right wing surface trailing edge up, radians 
BETA Angle of sideslip, positive nose right, radians 

P Roll rate, positive right wing down, radians/second 

PHI Bank angle, positive right wing down, radians 

R Yaw rate, positive nose right, radians/second 

RUD Rudder deflection, positive trailing edge Jeft, radians 


The SYSTEM POLES and DYNAMICS are listed the same as for page 3. 
Page 6 


This page lists all of the lateral-directional transfer functions in polynomial form for 
each state output to each contro! surface deflection input. It also gives the transfer 
function for the y-axis acceleration at the center of gravity to each control surface 
input. The variable names and units are the same as the state space model names 
and units on page 5. AY is the y-axis acceleration at the center of gravity. It is 
positive out the right side of the airplane and has units of feet/second?. 


References 


1. Heffley, R. K. and W. F. Jewell, "Aircraft Handling Qualities Data", NASA CR-2144, 
December 1972. 
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COMPUTER PRINTOUTS 


This section consists of the following Mach, altitude combinations: 


Mach Altitude Page 
(feet) 
3 100 13 
5 100 19 
7 100 25 
9 100 31 
95 100 37 
5 10000 43 
7 10000 49 
9 10000 55 
1.1 10000 61 
5 20000 67 
7 20000 73 
9 20000 79 
1.1 20000 85 
1.2 20000 91 
5 30000 97 
7 30000 103 
9 30000 109 
1.1 30000 115 
1.3 30000 121 
7 40000 127 
9 40000 133 
1.1 40000 139 
1.3 40000 145 
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